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Abstract. We investigate how different models that have been proposed for solving the dark matter problem can 
fit the velocity dispersion observed around elliptical galaxies, on either a small scale (~ 20kpc) with stellar tracers, 
such as planetary nebulae, or large scale (~ 200kpc) with satellite galaxies as tracers. Predictions of Newtonian 
gravity, either containing pure baryonic matter, or embedded in massive cold dark matter (CDM) haloes, are 
compared with predictions of the modified gravity of MOND. The standard CDM model has problems on a small 
scale, and the Newtonian pure baryonic model has difficulties on a large scale, while a fit with MOND is possible 
on both scales. 
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1. Introduction 

Measuring the velocity field in and around galaxies is the 
main way to test the dark matter distribution at small and 
intermediate scales. The observation of what are appar- 
ently non-Newtonian rotation curves around spiral galax- 
ies (e.g. Rubin et al. 1980) has been first solved by assum- 
ing that galaxies are embedded in dark matter haloes. 
Numerical simulations, however, predict a radial distribu- 
tion for the CDM model much more concentrated than 
what is observed (Gentile et al 2004; de Blok 2005). 

An alternative explanation was proposed by Milgrom 
(1983), as Modified Newtonian Dynamics (MOND). 
When the Newtonian acceleration falls below the critical 
value ao ~ 2 X 10^^" m s^^, the gravity law is empirically 
modified and then declines in l/r instead of 1/r^. Around 
spherical systems, the modified acceleration g satisfies the 
relation 

9^9 /aa) = 9n 

where gn is the Newtonian acceleration. For non-spherical 
geometry, this is only an approximation. However, in this 
Letter we only consider spherical systems as representing 
elliptical galaxies and we adopt ^{x) = x/Vl + x^. 

This model is very successful on a galactic scale; in 
particular, it explains a large number of rotation curves 
of galaxies, with some exceptions (Gentile et al. 2004), 
and naturally the TuUy-Fisher relation, (e.g. Sanders & 
McGaugh 2002). While the dark matter problem is ob- 
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servationally very clear around spiral galaxies, thanks to 
their rotation curve measured with the cold hydrogen gas 
at 21cm, which is nearly in circular orbits (e.g. Bosma 
1981, Verheijen & Sancisi 2001), the situation is much 
more complex around elliptical galaxies, with little or no 
rotation. 

Recently, planetary nebulae have been used as an ef- 
ficient tool for measuring the velocity field at large radii 
in early- type galaxies, and they complement stellar ab- 
sorption kinematical studies (Romanowky et al. 2003). In 
typical elliptical galaxies, the velocity-dispersion profiles 
were found to decline with radius, up to 5 effective radii, 
thereby requiring no dark matter at all. Dekel et al. (2005) 
show that the data are still compatible with the usual 
dark matter models, if the planetary nebulae tracers have 
particularly radial orbits in the outer parts, because of 
a recent merger with small impact parameter. However, 
the recent results from Douglas et al (2007) challenge this 
interpretation, prolonging the decline to more than 7 ef- 
fective radii. 

On larger scales around early-type galaxies, from 50 to 
SOOkpc, Klypin & Prada (2007), KP07, have proposed to 
test gravity models with satellite galaxies as a tracer. From 
the Sloan Digital Sky Survey, they stack several thou- 
sand galaxies in 3 luminosity classes and determine the 
number density of satellites and their velocity dispersion 
around them. In each mass range, the radial distributions 
are obtained with around 1500 satellites, although about 
1.5 satellite exist around each galaxy. 
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This large-scale galaxy neighborhood has not been 
widely tested yet in modified gravity. The well-known dif- 
ficulty of MOND in clusters has found a possible solution 
with neutrinos of 2 eV mass (Sanders 2003; Angus et al 
2007a) , and the escape velocity around giant galaxies like 
the Milky Way was shown to correspond to observations, 
when including the external field effect (Famaey et al 2007; 
Wu et al 2007). 

In this work, we solve the Jeans equation for the distri- 
bution of the velocity dispersion around elliptical galax- 
ies, and in particular we fit the NGC3379 galaxy, where 
the most extended data is available for the velocities. We 
also further explore fits at larger radii for the special case 
of NGC 3379, with satellite galaxies as tracers, as done 
by KP07 and Angus et al (2007b). This is only statisti- 
cally valid around a generic early-type galaxy, with a mass 
comparable to NGC 3379. This is the brightest galaxy of 
a group, but the observed companion velocity is not sta- 
tistically significant. 




10 20 

R (kpc) 



Fig. 1. Fit of the blue-light distribution in NGC 3379. 
The data points are taken from Douglas et al (2007). The 
curve represents the light profile derived with the mass- 
to-light ratio given in Table [T] and a mass model of cored- 
Hernquist type (Tabled]). 

2. The Jeans equation 

Because of the spherical symmetry, it has been shown 
(Angus et al. 2007) that the Jeans equation can be written 
in MOND, as well as in DM: 

— -I- (7 = -g{r) 

dr r 

where a is the radial velocity dispersion, a — dlnp/dlnr is 
the slope of the tracer density p, and /3 ~ 1 — (o'g+(T|)/2o'^ 
is the velocity anisotropy. 

The Jeans equation is usually used for a unique self- 
gravitating component, where the density p appearing in 
the "pressure" term on the left is the same as the density 
appearing in the Poisson equation, giving the density field 
g(r) on the right. On the left side, we use the density and 
the velocity dispersion of the tracers only, which can be 



Table 1. Parameters for deriving the mass of NGC 3379. 



D (Mpc) Ms Btc M/Lb 

9^8 ^109 10.10 ro~ 



Table 2. Radial distribution of the tracer density 
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very different from the density producing the potential, 
in particular for the satellite galaxies, which act as test 
particles. 

In the present approach, we want to compare all mod- 
els (Newtonian, with or without CDM and MOND) with 
the same dynamics for the Jeans equation, fitting the den- 
sity of tracers as close as possible to the observed density 
distribution. For the CDM models, we consider two differ- 
ent amounts of dark matter: (i) the CDMl model, reminis- 
cent of what is found in cosmological simulations, where 
the amount of dark matter inside the radius of 200kpc is 
equal to 60 times the visible matter; (ii) the CDM2 model, 
more akin to what is required to model rotation curves of 
spiral galaxies: the amount of dark matter inside the ra- 
dius of 200kpc is equal to 20 times the visible matter. 

All models are then compared with the same tools, 
while KP07 compare the CDM model with numerical 
simulations, and MOND with the Jeans equation, but 
without exploring all variations of the tracer density and 
anisotropy profiles. 

3. Results 

NGC 3379 is modeled by a cored-Hernquist distribution 
(see Fig. [T|). The stellar mass of this galaxy is estimated 
to be M = 1 X lO^^AfQ, assuming a constant mass-to-light 
ratio M/L (see Table [1]). In MOND (and Newton without 
DM) , this stellar distribution is the only matter contribut- 
ing to the gravitational potential. In Newtonian dynamics 
with DM, the DM halo follows the NFW radial profile 
found in N-body cosmological simulations. We used three 
kinds of tracers to compare Newton without DM, CDMl, 
CDM2, and MOND models, for different length scales: 
stars for the inner part (r < 1 kpc, Shapiro et al 2006), 
planetary nebulae for the middle part (1 kpc < r < 20 
kpc, Douglas et al. 2007), and satelhtes (20 kpc < r < 200 
kpc, KP07). The satellites do not correspond to the real 
ones orbiting NGC 3379, but are a statistical representa- 
tion from the SDSS (KP07), giving the velocity dispersion 
of satellites submitted to the gravity of a typical ellipti- 
cal galaxy in the same mass category as NGC 3379. The 
parameters of the tracer density are displayed in Table [2l 
We numerically solve the Jeans equation for each 
tracer. The projected velocity dispersion (ctjos) is repre- 
sented in Fig. [21 
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Fig. 2. Fit of the velocity dispersion data around NGC 3379 with pure a Newton model without DM, CDMl, CDM2, 
and MOND. The data points in the very center (r < 1 kpc) are from Shapiro et al. (2006). On intermediate scales 
they are taken from Douglas et al. (2007) and are based around 10-30kpc on the planetary nebulae as tracers. On 
large scales, the data points correspond to the statistical distribution from the satellite galaxies as tracers (KP07), 
adapting NGC 3379 in their small luminosity class. The tracer density is taken as an Hernquist distribution on a 
small scale, and fitted to the satellite density distribution on a large-scale. The best fit is obtained with varying /3(r). 
Purple/dot-dash: CDMl. Green/dash: CDM2. Blue/fuU line: MOND, Cyan/short dash: Newton without DM. 



3.1. Newton dynamics without DM 

— Stars: the projected velocity dispersion ((T*^^) is well- 
fitted using a small radial anisotropy increasing with 
radius from f3 ~ 0.2 to /3 = 0.4. Since the acceleration 
here is above the critical oq, both MOND and Newton 
share the same fitting parameters. 

— FN: as for the stars, the observed FN velocity dis- 
persion {(jfj^ ) is also in good agreement with the 
Newtonian model. It implies a tangential velocity in 
the central region {(3 = —0.6) evolving to radial tra- 
jectories {P = 0.8). 

— Satellites: af^l cannot be fitted at all by considering 
the stellar mass alone. A dark matter component is 
required. 



3.2. CDMl: DM halo mass from cosmological 
simulations 

The mass of the NFW halo is 60 times the stellar mass 
inside 200fcpc. 



— Stars: the cusp of the DM halo does not allow fit- 
ting the (7*^^, whatever the anisotropy profile. The pre- 
dicted velocity dispersion is too large. 

— FN: (Tjq^ does not fit the entire data set. The disper- 
sion is too large around 10 kpc. If (3 is increased to- 
wards its maximum, up to a purely radial anisotropy 
profile (which tends to decrease the velocity disper- 
sion), af^^ is then too small in the inner part. In the 
best fit, the value of (3 is maximum at 1 kpc. 

— Satellites: the model reproduces the observations well, 
if radial anisotropy is maintained all over the region 
iP > 0.8). 

3.3. CDM2: DM halo mass compatible with average 
spirals 

The mass of the NFW halo is 20 times the stellar mass 
inside 200fcpc. 

— Stars & FN: they are well-fitted, and the DM halo 
cusp is no longer a problem in the internal dynamics 
of the galaxy. For the FN, the velocity dispersion is 
isotropic near the center (r < 1 kpc) and increase to 
its maximum /3 = 1 at 20 kpc. 
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— Satellites: it is more difficult to fit crf^^ with less DM. 
We need to impose an isotropic velocity dispersion 
(/3 ^ 0), to avoid too fast Oa fall in af^^. The nega- 
tive gradient of /3{r) helps to straighten the slope of 
(Ji^Kr) which otherwise is too steep. 

We also varied the dark-to-visible mass ratio between 
10 to 60, and 20 gives the best compromise between the 
1 — 10 kpc and 30 — 200 kpc regions. 

3.4. MOND 

All scales are easy to fit with MOND. 

— Stars: under 1 kpc, the gravitational potential is purely 
Newtonian, so MOND and Newton models are iden- 
tical. The anisotropy needs to be a bit radial {(3 ~ 
0.2 — 0.4), but always far from the allowed limits. 

— PN: until 5 kpc the potential is still Newtonian. A 
slight difference appears after 8 kpc. When varying 
/3(r) from -0.6 (tangential) to 0.8 (radial), the MOND 
model is in good agreement with the observations. 

— SateUites: a^^l is also well-fitted by MOND on this 
scale, using /3 = —0.2 at r = 20 kpc to /3 = 0.2 at 
r = 200 kpc. 

Figure [3] shows how the variations in (3 can help to 
fit many versions of the velocity curves. As in all other 
models, there is still much latitude to the fit. Large spatial 
variations of /3 are possible if the ellipticals are the results 
of mergers (Dekel et al., 2005) 

4. Discussion 

In contrast to the conclusion of KP07, MOND does not 
predict constant velocity dispersion with radius in the 
neighborhood of elliptical galaxies. There is a wide lat- 
itude for varying the possible anisotropy parameter ac- 
cording to the scale and the tracer considered. Our best 
fit starts from a tangential or isotropic configuration near 
the center and evolves progressively to a radial one for 
each of the three scales considered, which appears quite 
realistic. 

While the Newtonian model without dark matter has 
problems in the outer parts, the CDMl model encounters 
severe difficulties in the inner parts. The CDM2 model. 
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Fig. 3. Fit of the velocity dispersion data around NGC 
3379 with MOND, and various values of j3. The data 
points are the same as in Fig. [51 The 4 curves on a small- 
scale are 0., 0.5, 0.9; and (3=-l, 0.,0.5 on a large- 
scale. 

with a reduced dark matter relative to the visible mass, 
can also fit the data quite well. Compared to the MOND 
model, it requires a larger radial anisotropy on each scale, 
and its /3 profile is unusual for the satellite tracer with a 
negative slope. For /3= constant, the crf°^* would be too 
steep compared to the observations. And by introducing 
anisotropy increasing with radius, the (Ti^^ slope is even 
more increased. 

In the MOND regime, the external field effect (EFE, 
e.g. Wu et al 07) would also change the predictions. This 
will modulate the actual force on the particle tracers, and 
on the velocity. Since there was still latitude for fitting the 
observations with MOND, we feel that it is still possible 
to consider it with other anisotropy parameters. However, 
we note that the external field effect is not even necessary, 
for reproducing the observed velocity dispersion slope. 
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